Lorentz covariant nucleon self-energy decomposition of the nuclear symmetry energy 
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Using the Hugenholtz-Van Hove theorem, we derive analytical expressions for the nuclear symme- 
try energy E sym (p) and its density slope L(p) in terms of the Lorentz covariant nucleon self-energies 
in isospin asymmetric nuclear matter. These general expressions are useful for determining the 
density dependence of the symmetry energy and understanding the Lorentz structure and the mi- 
croscopic origin of the symmetry energy in relativistic covariant formulism. As an example, we 
analyze the Lorentz covariant nucleon self-energy decomposition of E sym (p) and L(p) and derive the 
corresponding analytical expressions within the nonlinear a-u-p-5 relativistic mean field model. 
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I. INTRODUCTION 



In the current research of nuclear physics and astro- 
physics, there is of great interest to study the density de- 
pendence of the nuclear symmetry energy E sym (p) that 
essentially characterizes the isospin dependent part of 
the equation of state (EOS) of asymmetric nuclear mat- 
ter. The exact knowledge on the symmetry energy is 
important for understanding not only many problems in 
nuclear physics, but also many critical topics in astro- 
physics pHll as well as some interesting issues regarding 
possible new physics beyond the standard model 0-0- 
During the last decade, although significant progress has 
been made both experimentally and theoretically on con- 
straining the density dependence of the symmetry en- 
ergy @, 1, EMI ( see > c -g-> Refs - EM3 for review 
of recent progress), large uncertainties on E sym (p) still 
exist, especially its super- normal density behavior re- 
mains elusive and largely controversial [19l - |22l ]. To re- 
duce the uncertainties of the constraints on E sym (p) is 
thus of critical importance and remains a big challenge 
in the community, and this provides a strong motivation 
for studying isospin nuclear physics in radioactive nuclei 
at the new/planning rare isotope beam facilities around 
the world, such as CSR/Lanzhou and BRIF-II/Bcijing 
in China, RIBF/RIKEN in Japan, SPIRAL2/GANIL in 
France, FAIR/GSI in Germany, FRIB/NSCL in USA, 
SPES/LNL in Italy, and KoRIA in Korea. 

Theoretically, the uncertainties of the constraints on 
E sym (p) are mainly due to the lack of knowledge about 
the isospin dependence of in-medium nuclear effective in- 
teractions and the limitations in the techniques for solv- 
ing the nuclear many body problem. Very recently, it has 
been proposed that it is very useful to directly decompose 
E sy m(p) in terms of some relevant parts of the commonly 
used underl ying nuclear effective interaction in nuclear 
medium 0M[24|. Based on the Hugenholtz-Van Hove 
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(HVH) theorem, indeed, the E sym (p) can be decomposed 
analytically in terms of the single-nucleon potential in 
asymmetric nuclear matter and the resulting expressions 
are quite general and independent of the detailed nature 
of the nucleon interactions, providing an important and 
physically more transparent approach to extract informa- 
tion on the symmetry energy from the isospin dependence 
of strong interaction in nuclear medium and understand 
why the symmetry energy predicted from various models 
is so uncertain [22|,[24[. In these works, the decomposi- 
tion of -E S ym (p) is based on non- relativistic framework. It 
is thus of great interest to explore more general decom- 
position within relativistic covariant framework, which is 
the main motivation of the present work. 



II. COVARIANT SELF-ENERGY 
DECOMPOSITION OF E Bym (p) AND L{p) 

The relativistic covariant formulation has made great 
success during the last decades in understanding many 
nuclear phenomena [25l - l27| . In particular, the mi- 
croscopic relativistic covariant Dirac-Brucckncr-Hartree- 
Fock (DBHF) approach (28rl33j has achieved impressive 
success in describing the saturation properties of nuclear 
matter without any need to introduce a three-nucleon 
force required in the microscopic non-relativistic BHF 
calculations (see, e.g., Refs. [HI, HH). It has been argued 
that in non-relativistic calculations the three-nucleon 
forces must be introduced to mimic the variation of the 
Dirac spinors in the nuclear medium contained in rela- 
tivistic covariant approach (32j . In addition, the Lorentz 
covariant decomposition of the nuclear mean field poten- 
tial has been shown to be very important for understand- 
ing the dynamics in heavy ion collisions at relativistic en- 
ergies |36H38j . These features imply that the Lorentz co- 
variance could be important for understanding the higher 
energy /density nuclear phenomena, e.g., the high density 
behaviors of the symmetry energy. 

Owing to the translational, rotational and time- 
reversal invariance, parity conservation, and hcrmitic- 
ity, the Lorentz covariant nucleon self-energy in the rest 
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frame of asymmetric nuclear matter with baryon density 
P = Pn + Pp and isospin asy mmetry a = (p n — p p )/p can 
be written generally as [25], [H, [39M4lT | , 

£ J (p, a, |k|) =^(p, a, |k|) - y^> J (p, a, |k|) 



=£^,a,|k|)+ 7 °5^(p,a,|k|) 



(1) 



+ 7 -k°S 7 (p,a,|k|), 

where J = n or p is isospin index; Sg(/0, a, |k|) is 
the scalar self-energy, Ey(p,a, |k|) = — TP' J {p, a, |k|) 
and Sj^(/0, a, |k|) are, respectively, the zeroth compo- 
nent and the space component of the vector self-energy 
E A '' J (/9, a, |k|), and they all generally depend on p, a and 
the magnitude of the nucleon momentum |k| (Minkowski 
metric is = (+, — , — , — )); k° = k/|k| is the unit vec- 
tor of momentum k. A proof of Eq. ([T]) can be found in 
Ref. 25 1. Accor dingly, the single-nucleon energy can be 
expressed as [UHlMEl 



£ J (p, a, |k|) = £* J (p, a, |k|) + £^(p, a, |k|), 



(2) 



with 



£* J {p, a, |k|) = yj\k*'(p, a, \k\)Y + [M}(p, a, \k\)Y, 

(3) 

where the nucleon effective (Dirac) mass Mj(p, a, |k|) 
and effective momentum k J *(p, a, |k|) are defined, re- 
spectively, as 



A/}(p,a,|k|) = M + Hi(p,a,\k\), 
k* J (p,a,|k|) = k + k°S^(p,a,|k|), 



(4) 
(5) 



with M being the nucleon mass. 

Due to the exchange symmetry between protons and 
neutrons in nuclear matter, the EOS of asymmetric nu- 
clear matter, defined by its binding energy per nucleon, 
can be expanded as a power series of even-order terms in 
a as 



E{p,a) ~ E (p) + E sym (p)a 2 + 0(a 4 ), 



(6) 



where Eq( P ) = E(p,a = 0) is the EOS of symmetric 
nuclear matter, and the symmetry energy is expressed as 



E. 



sym(p) 2 



1 d 2 E{p,a) 



da 2 



(7) 



a=0 



Around the nuclear matter saturation density pa, the 
symmetry energy can be expanded as 



E syia (p) ~ E sym (p ) + L X + 0{ X 2 ), 



(8) 



where x = (p ~ Pa)/ 5 Pa is a dimcnsionlcss variable, and 
L = L(po) is the density slope parameter of the symme- 
try energy at po and thus carries important information 
on the symmetry energy at both high and low densities. 
Generally, the slope parameter of the symmetry energy 
at arbitrary density p is defined as 



HP) 



3p ^^ s y m 

dp 



(p) 



(9) 



According to the HVH theorem [42|, |43| , the nucleon 
chemical potential in asymmetric nuclear matter should 
be equal to its Fermi energy (the single-nucleon energy 
at Fermi surface), i.e., 



£^{p, a, fcfO 



d[pE{ Pl a)} 
dpj 



M, 



(10) 



where £p(p, a,fcp) = £ J (p,a,kp) is the nucleon Fermi 
energy, and fcp = kp(l + r/a) 1 / 3 (we assume t£ = 1 
and Tg = — 1 in this work) is the nucleon Fermi momen- 
tum with kp = (37r 2 p/2) 1 / 3 being the Fermi momentum 
in symmetric nuclear matter at density p. It should be 
noted that the HVH theorem is independent of the de- 
tailed nature of the interactions used and is valid for any 
interacting self-bound infinite Fermi system, such as in- 
finite nuclear matter (4l| - |43j . 

Expanding E(p, a) as a power series of even-order 
terms in a on the right-hand side of Eq. (TIT))) , we can 
obtain 

ri£i(p,a,k^) = Y^^E sym ,2 l (p)a 2l -\ (11) 

J— p,n i— 1 

£ £i( P ,a,ki)=2^^ + 2M 



J—p,n 



dp 

+^ p Y.p E ^^p) a21 ^ ( 12 ) 



where E sym , 2i (p) = [l/(2i)\}[d 2z E{p,a)/da 2 % =0 are the 
symmetry energies of different orders and particularly 
we have E SYm ^{p) = E sym (p). Furthermore, expanding 
£p(p,a, fcp) as a power series of a on the left-hand side 
of Eq. (fTT|) and Eq. (|T2"j) , and comparing the coefficients 
of the first-order a terms on both left- and right-hand 
sides of Eq. ([TTj) , we then obtain 



Esym(p) 



1 d 
4 da 



J2 -riH{p,<x,H) 

J=p,n 



(13) 



Q=0 



while comparing the coefficients of second-order a terms 
on both sides of Eq. (fT2j) leads to the following expres- 
sion: 



3d 2 
L{p) = Ada>- 



E 

J— p,n 



£p(p,a,kp) 



+ 3E S 



a=0 



^sym(p)- 

(14) 

Substituting Eq. © into Eq. (|T3)) . we can obtain 



sym 

_ zriOjmom.K 
sym 

_ 77ilst,K 
sym 



(p) 



^0,mom,S / 
;ym 

! (p) 



£ s u " 5 (p) 



0,mom,V 



E sym VP) + fi sym 



E' 

sym 
(Pi 



(P) 

(15) 



where E^(P), E^^°(p) and E^f(p) (here O de- 
notes K, S or V) represent, respectively, the contribu- 
tions from the kinetic part, the momentum dependence 
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of the nucleon self-energies in symmetric nuclear matter 
and the first-order symmetry self-energies, and they can 
be expressed analytically as 



£ s °" K (p) 

(p) 
(p) 

(p) 
(p) 



k F k* F ( P ) 0E°(p,|k|) 



6^( P ) a|k| 

>M *(p) 0Eg(p,|k| 



E, 



syr 

0,mom,V / 
sym 

irilst.K/ 

sym 



6£|(p) d|k 
_fc F 9S°(p,|k|) 



|k|=fe F 



k| — fcp 



6 



d|k| 



k|=/SF 



"2 £ F (p) 
l M *(p)£ s / m4 (p,|k| 

"2 
1. 



fc F ) 



sym 



^sym,l 



£p(p) 

(p, |k| = fc F ), 



(16) 
(17) 
(18) 
(19) 
(20) 

(21) 
(22) 



fc F + ££(p, fc F ), M *(p) = M + £°(p, fc F ), 
-M* 2 )V2, £»(p,|k|) = ££(p,a = 0,|k|), 
E£(p,a = 0,|k|), E^(p,|k|) = = 
0, |k|), and the «-th order symmetry self-energy is defined 
as (here O = K,S,V) 



where fc F (p) = 

£f(p) = (k? 
Sg(p,|k|) 



^ T 3 J 'SJ(p,q,|k|) 



J— p,n 



Q=0 

(23) 



Furthermore, Eq. f| 15|) can be rewritten as 



-Esym (p) 



6M, 



O.Lan 



(P,|k 



E, 



lst,K 
sym 



|k|=fc 
lst.S 



(p) + E^(p) + E^(p), (24) 



where Afg Lan (p, |k|) is the nucleon Landau mass 
in symmetric nuclear matter, i.e., Afg Lan (p, |k|) = 
\k\[d\k\/d£°(p, IkQ]- 1 (see, e.g., Ref.' 0) with 
£°(p, |k|) = £ J (p,a = 0, |k|), and one can 
easily verify the relation fc F /6Mg* Lan (p, fc F ) = 

£^»(p) + £ s ° y £° m ' K (p) + Ss° y S° m ' S (p) + £s° ym ° m ' V (p)- In 
this way, we have decomposed analytically the symmetry 
energy E sym (p) in terms of the Lorentz covariant nucleon 
self-energies in asymmetric nuclear matter. 

Similarly, by substituting Eq. ([2]) into Eq. (fi4|) . the 
slope parameter L(p) can be decomposed as 



L(P) 



=L kin (p) + L mom (p) + L lst (p) 
+ L cross (p) + L 2nd (p), 



with 



L kin (p) 



k F k F 
~6£T 



6£f ' 



(25) 



(26) 



\P) 



/,- 2 

F 



3£ F 2 

fcpa 2 ^ 



d|k| 



Mq # 2 e° 



a 2 s° 



fc 2 , 

67^ 



£p <9|k| 2 £* d\k\ 2 d\k\ 



d\k\ 



Ikl — Aif 



£p S|k 



£ F 5|k| 



iki=fe 

2 



J |k|=fcp 



fcp fcp Afg 
3£f 



d|k| 



d|k| 



(27) 



L lst (p) 



2£*3 



*yisym,l 



M n *E 







fcp£ s 



k|=fe F 
sym,l 



3 
2 



F v^ s y m -i 

K 



y^sym.l 
c* S 
C F 



risym. 1 



2 v sym,l 



fcp/cpAfg £ g y 111,1 



c*3 
C F 



\P) = 



k* as^" 1 ' 1 



£ F 0|k 


i 






°F 


£f \ 


fc F £ s s ym ' 1 


rM* 2 




c*2 
L C F 



<?*3 
C F 



(28) 



9|k| 



K d\k\ 



k F <9£k 



L 2nd (p) 



M * 9|k| 0|k| 



^F v sym,2 . j^Q v sym,2 
~c* K "+■ T^T^S 
^•p ^p 



J lkl= 




9£°1 




0|k|J 


|k|=fep 


9£0 




9|k| 


ik|=fe P 



v^sym.2 



(29) 
(30) 



On the right-hand side of Eqs. (|2"6 > (|3"0 ]) . the density 
and momentum dependence have been suppressed with 
^ m,i (p, |k| = fc F ) (O = K,S,V). Eq. (Jig (or 
and Eq. (|2"5"|) are two main results of this work. 



v^sym,i 



III. APPLICATION TO THE NONLINEAR 
cT-u-p-5 RMF MODEL 

The nucleon self-energies can be calculated theoreti- 
cally from a certain relativistic covariant approach or ex- 
tracted experimentally (around po) from the Dirac phe- 
nomenology of nucleon-nucleus scattering. The Lorentz 
covariant nucleon self-energy decompositions of E syTCL (p) 
in Eq. ([ISJ (or ([21])) and L{p) in Eq. (gSJ) are gen- 
eral and they are useful for determining the density de- 
pendence of the symmetry energy and understanding its 
Lorentz structure and the microscopic origin. As an ex- 
ample, we consider here the nonlinear a-uj-p-8 relativistic 
mean field (RMF) model which is based on effective inter- 
action Lagrangians involving nucleon and meson fields, 
and has been widely discussed in the literature (see, e.g., 
Ref. [3]). A very useful feature of this model is that 
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the nucleon self-energies in asymmetric nuclear matter 
can be obtained analytically and this makes our analysis 
physically transparent. The Lagrangian density of the 
nonlinear a-uj-p-5 RMF model can be expressed as (see, 
e.g., Ref. 0): 

C =i> fr^ifl" - gu u») - (M - g a a)\ $ 



1 



1 



1 



b a M( g<7 a) 3 - 7 c a (g a a) 4 + -c^a^) 



p,\2 



+ o (9 2 P Pn ■ P M ) (A s .9> 2 + AvflSw/iW) 



9 P Pp ■ ^f J i) + g s 8 ■ ipr J i), 



(31) 



where F pu = d^ v - d„u p and G M „ = d il p v - d v p^ 
are strength tensors of uj field and p field, respectively 
while V; Pp an d S arc nucleon field, isoscalar- 

scalar field, isoscalar-vector field, isovector-vector and 
isovector-scalar field, respectively, and the arrows denote 
the isospin vector. As and Ay are two cross-coupling 
constants for varying the density dependence of E sy m(p), 
and m CT , m u , m p , and m$ are masses of mesons. 

In the RMF model, meson fields arc replaced by their 

— (3) -* 

expectation values, i.e., a — > a, cuq — > p — > /9„, 
where the subscript "0" denotes the zeroth component 
of the four-vector while the superscript "(3)" denotes the 
third component of isospin, Furthermore, the space-like 
self-energy 8, |k|) vanishes (due to the Hartree ap- 

proximation in the RMF model) while the scalar and 
time- like self-energies are momentum independent, i.e., 



(p,a) 



£ v (/>,a) 



6<- 3 \ 



-g a a + t 3 g s 



J -(3) 

g^ujQ ~ t 3 g p pl . 



(32) 
(33) 



The symmetry energy then can be decomposed as 

pkin / \ , z7'lst,S/„\ i 7^1st,V/ 



E sym (p) 



zE sym(p) + Kymip) + E sym (p) 



2 si 



1. 



^sym,l 



(p), (34) 



where the (lst-order) symmetry self-energies are 



^sym,l 



isym,l 



(P) 



(p) 



9 2 P P 

Qp 1 



(35) 
(36) 



with Q s = m 2 + 3g$(ps/MQ — p/£p), ps being the scalar 



density, and Q p 



+ ^sgigp 2 + Avgigivi. We 



note that the above analytical expression for E sym (p) is 
exactly the same as the one obtained from the normal 
approach (see, e.g., Ref. [H[)- Similarly, the slope pa- 
rameter L(p) can be decomposed as 



L(p) = L^(p) + L^(p)+L^(p) 



(37) 



with 



L kin (p) 



k 2 (£* 2 +Mf 



L lst ( P ) = r 



6£* 3 
M *E° yDM 



(38) 



(p) 



r-isym.,1 



(p) 



3fcp 



2£* 3 



ipsym.l 



L 2nd (p) =3 



A/ *E s s ym ' 2 



(p) 
(P) 



2 M *^E^ m ' 1 



(p) 



(7*3 



+ S^ m ' 2 (p) 



, (39) 
(40) 



where the 2nd-order symmetry self-energies are 



ysym,2 



(p) 



7 |M ' a P r 



2^A/ ^|p 2 



2Q 



£p 2 Qs 



v»sym,2 



(p) 



g a M*k 2 p 
3£* 3 

9^9 4 p A ^oP 2 



£p 4 Qs 

2glg 4 p A s ap 2 \ 



(41) 



(42) 



m' 



gl(3 Ps /MZ ~ Z P /£* F ) + IKMgla - 
^ Zc^giul and T = 3 5(t5 2 (2 Ps /M * 2 - 



with Q CT = 
3c CT ^a- 2 ,<; 

3p/M££ F + M£p/£f). To the best of our knowledge, 
the above formulas give, for the first time, the analyti- 
cal expression of the slope parameter L(p) in the non- 
linear RMF model. The above analytical expressions of 
E sym {p) and L{p) can be easily generalized to the case 
of the density dependent RMF model that has similar 
isospin structure as the nonlinear RMF model 44(. It 



should be mentioned that these analytical expressions 
for E sym (p) and L(p) are very useful for determining the 
isovector parameters in the RMF model by fitting the 
empirical prop erties of asymmetric nuclear matter (see, 
e.g., Ref. j45[ for such a procedure in the case of the 
isoscalar sector). 

Shown in Fig. [T]is the density dependence of E sym (p) 
and its self-energy decomposition according to Eq. 
for four interactions^ i.e., FSUGold 0, IU-FSU gj, 
NLpc5 H3 and HA gg|. FSUGold and IU-FSU are two ac- 
curately calibrated interactions based on the ground state 
properties of closed-shell nuclei, their linear response, and 
the structure of neutron stars. Since FSUGold and IU- 
FSU do not consider the isovector-scalar 6 meson field, 
one thus has Eg ym,1 (p) = 0. On the other hand, both 
NLpS and HA include the S meson field with the latter 
fitting successfully some results obtained from the micro- 
scopic DBHF approach while the former fitting the em- 
pirical properties of asymmetric nuclear matter and de- 
scribing reasonably well the binding energies and charge 
radii of a large number of nuclei [H[ . From Fig. [TJ one can 
see that while the kinetic contribution E^f n (p) is roughly 
the same for different interactions, the different interac- 
tions predict significantly different values for E^^(p) 
and Egy^(p). In particular, compared with FSUGold 
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FIG. 1: (Color online) Density dependence of E ayul (p) and its 
self-energy decomposition according to Eq. (|34|) in the non- 
linear RMF model with different interactions. 
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FIG. 2: (Color online) Density dependence of L(p) and its self- 
energy decomposition according to Eq. (|37[) in the nonlinear 
RMF model with different interactions. 



and IU-FSU, HA predicts very similar total E sym (p) but 
significantly different Ej^(p) and E^^(p). 

Similarly, we show in Fig. [2] the density dependence 
of the slope parameter L(p) and its self-energy decom- 
position according to Eq. (JSTJ) for FSUGold, IU-FSU, 
NLpS and HA. Again, it is seen that the different inter- 
actions predict roughly same kinetic contribution L kln (p) 
but significantly different values for L lst (p) and L 2nd (p). 
In particular, one can see that the higher-order con- 
tribution L 2ad (p) from the second-order symmetry self- 
energies generally cannot be neglected, agreeing well with 
the recent non-relativistic calculations I24II. 



IV. SUMMARY AND OUTLOOK 

Using the Hugenholtz-Van Hove theorem, we have 
shown that the symmetry energy E sym (p) and its den- 
sity slope L{p) can be decomposed analytically in terms 
of the Lorentz covariant nucleon self-energies in asym- 
metric nuclear matter, and the corresponding expressions 
have been derived for the first time. These general ex- 
pressions for the covariant self-energy decomposition of 
E sym {p) and L(p) are useful for determining the den- 
sity dependence of the symmetry energy, deciphering the 
Lorentz structure of the symmetry energy, and under- 
standing the microscopic origins of the symmetry energy. 
As an example, we have analyzed the Lorentz covariant 
nucleon self-energy decomposition of E sym (p) and L{p) 
within the nonlinear a-tu-p-S relativistic mean field model 
and derived the corresponding analytical expressions for 
E sym (p) and L(p), which are potentially useful for fixing 
the isovector parameters in the RMF model from fitting 
the empirical properties of asymmetric nuclear matter. 



From analyzing the self-energy decomposition of 
-E'sym(p) and L(p) within the nonlinear a-uj-p-S rela- 
tivistic mean field model, we have found that the re- 
sults strongly depend on the interactions used and also 
whether the isovcctor-scalar S meson is included or not. 
These results imply that it is of great importance to de- 
termine individually each part of the Lorentz covariant 
nucleon self-energy decomposition of E sym (p) and L(p) 
from experiments (e.g., Dirac phenomenology ) or mi- 
croscopic calculations based on nucleon-nucleon interac- 
tions derived from scattering phase shifts (e.g., DBHF). 
On the other hand, the Lorentz covariant nucleon self- 
energies in asymmetric nuclear matter can also be deter- 
mined from quantum chromodynamics (QCD) by means 
of QCD sum-rule techniques [50Ll5ll ]. The general expres- 
sions of the Lorentz covariant nucleon self-energy decom- 
position of -E'sym(p) and L(p) presented in this work are 
thus very useful for determining the symmetry energy 
from QCD. These studies are in progress. 

Acknowledgments 

The authors would like to thank R. Chen, B.A. Li, 
X.H. Li and C. Xu for useful discussions. This work was 
supported in part by the NNSF of China under Grant 
Nos. 10975097 and 11135011, the Shanghai Rising- 
Star Program under Grant No. 11QH1401100, the "Shu 
Guang" project supported by Shanghai Municipal Ed- 
ucation Commission and Shanghai Education Develop- 
ment Foundation, the Program for Professor of Special 
Appointment (Eastern Scholar) at Shanghai Institutions 
of Higher Learning, the Science and Technology Com- 
mission of Shanghai Municipality (11DZ2260700), and 



6 



the National Basic Research Program of China (973 Pro- gram) under Contract No. 2007CB815004. 



B.A. Li, CM. Ko, and W. Bauer, Int. Jour. Mod. Phys. 
E 7, 147 (1998). [25 
P. Danielewicz, R. Lacey, and W.G. Lynch, Science 298, 
1592 (2002). [26 
J.M. Lattimer and M. Prakash, Science 304, 536 (2004); 
Phys. Rep. 442, 109 (2007). [27 

A. W. Steiner, M. Prakash, J.M. Lattimer, and P.J. Ellis, 
Phys. Rep. 411, 325 (2005). [28 
V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. 
Rep. 410, 335 (2005). [29 

B. A. Li, L.W. Chen, and CM. Ko, Phys. Rep. 464, 113 
(2008). [30 

C. J. Horowitz, S.J. Pollock, P.A. Souder, and R. 
Michaels, Phys. Rev. C 63, 025501 (2001). [31 
T. Sil, M. Centelles, X. Vinas, and J. Piekarewicz, Phys. 
Rev. C 71, 045502 (2005). [32 

D. H. Wen, B.A. Li, and L.W. Chen, Phys. Rev. Lett. 

103, 211102 (2009). [33 
L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. Lett. 94, [34 
032701 (2005); Phys. Rev. C 72, 064309 (2005); B.A. Li 
and L.W. Chen, Phys. Rev. C 72, 064611 (2005). [35 
M.B. Tsang, Y. Zhang, P. Danielewicz, M. Famiano, Z. 
Li, W. G. Lynch, and A. W. Steiner, Phys. Rev. Lett. [36 
102, 122701 (2009). 

M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, [37 
Phys. Rev. Lett 102, 122502 (2009); M. Warda, X. 
Vinas, X. Roca-Maza, and M. Centelles, Phys. Rev. C [38 
80, 024316 (2009). 

J.B. Natowitz, G. Ropke, S. Typel, D. Blaschke, A. [39 
Bonasera, K. Hagel, T. Klahn, S. Kowalski, L. Qin, S. 
Shlomo, R. Wada, and H. H. Wolter, Phys. Rev. Lett [40 

104, 202501 (2010). 

C. Xu, B.A. Li, and L.W. Chen, Phys. Rev. C 82, 054607 [41 

(2010) . 

L.W. Chen, CM. Ko, B.A. Li, and J. Xu, Phys. Rev. C [42 
82, 024321 (2010). 

M.B. Tsang, Z. Chajecki, D. Coupland, P. Danielewicz, [43 
F. Famiano, R. Hodges, M. Kilburn, F. Lu, W.G. Lynch, 
J. Winkelbauer, M. Youngs, and Y.X. Zhang, Prog. Part. [44 
Nucl. Phys. 66, 400 (2011). 

L.W. Chen, Phys. Rev. C 83, 044308 (2011). [45 
W.G. Newton, M. Gearheart, J. Hooker, and B.A. Li, 
larXiv:1112.20T8l [46 
Z.G. Xiao, B.A. Li, L.W. Chen, G.C Yong, and M. 
Zhang, Phys. Rev. Lett. 102, 062502 (2009). [47 
Z.Q. Feng and CM. Jin, Phys. Lett. B683, 140 (2010). 
P. Russotto, P.Z. Wu, M. Zoric, M. Chartier, Y. [ 
Leifels, R.C Lemmon, Q. Li, J. Lukasik, A. Pagano, P. 
Pawlowski, and W. Trautmann, Phys. Lett. B697, 471 [49 

(2011) . 

C. Xu and B.A. Li, Phys. Rev. C 81, 064612 (2010). [50 
C. Xu, B.A. Li, L.W. Chen, and CM. Ko, Nucl. Phys. 
A865, 1 (2011). [51 
[24] R. Chen, B.J. Cai, L.W. Chen, B.A. Li, X.H. Li, and C. 



Xu, Phys. Rev. C 85, 024305 (2012). 

B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 

(1986) . 

B. D. Serot and J.D. Walecka, Int. J. Mod. Phys. E 6, 515 
(1997). 

J. Meng, H. Toki, S.G. Zhou, S.Q. Zhang, W.H. Long, 
and L.S. Ceng, Prog. Part. Nucl. Phys. 57, 470 (2006). 
MR. Anastasio, L.S. Celenza, W.S. Pong, and CM. 
Shakin, Phys. Rep. 100, 327 (1983). 

C. J. Horowitz and B.D. Serot, Nucl. Phys. A464, 613 

(1987) . 

R. Brockmann and R. Machleidt, Phys. Rev. C 42, 1965 
(1990). 

E.N.E. van Dalen, C. Fuchs, and A. Faessler, Eur. Phys. 
J. A 31, 29 (2007). 

E. N.E. van Dalen and H. Miither, Int. J. Mod. Phys. E 
19, 2077 (2010). 

F. Sammarruca, Int. J. Mod. Phys. E 19, 1259 (2010). 
Z.H. Li, U. Lombardo, H.-J. Schulze, W. Zuo, L.W. Chen, 
and H.R. Ma, Phys. Rev. C 74, 047304 (2006). 

Z.H. Li, U. Lombardo, H.-J. Schulze, and W. Zuo, Phys. 
Rev. C 77, 034316 (2008). 

CM. Ko, Q. Li, and R. Wang, Phys. Rev. Lett. 59, 1084 
(1987). 

B. Blattel, V. Koch, and U. Mosel, Rep. Prog. Phys. 56, 
1 (1993). 

T. Gaitanos, M. Di Toro, S. Typel, V. Baran, C. Fuchs, 
V. Greco, and H.H. Wolter, Nucl. Phys. A732, 24 (2004). 
M. Jaminon, C. Mahaux, and P. Rochus, Nucl. Phys. 
A365, 371 (1981). 

C. J. Horowitz and B.D. Serot, Nucl. Phys. A399, 529 
(1983). 

H. Uechi, Phys. Rev. C 41, 744 (1990); Nucl. Phys. 
A696, 511 (2001). 

N.M. Hugenholtz and L. Van Hove, Physica 24, 363 
(1958). 

L. Satpathy, V.S. Uma Maheswari, and R.C. Nayak, 
Phys. Rep. 319, 85 (1999). 

L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. C 76, 
054316 (2007). 

N.K. Glendenning, Compact Stars, 2nd Edition, 2000, 
Spinger-Verlag New Yrok, Inc., p. 178. 
B.G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett. 
95, 122501 (2005). 

F.J. Fattoyev, C.J. Horowitz, J. Piekarewicz, and G. 
Shen, Phys. Rev. C 82, 055803 (2010). 
B. Liu, V. Greco, V. Baran, M. Colonna, and M. Di Toro, 
Phys. Rev. C 65, 045201 (2002). 

J.K. Bunta and S. Gmuca, Phys. Rev. C 68, 054318 
(2003). 

T.D. Cohen, R.J. Furnstahl, and D.K. Griegel, Phys. 
Rev. Lett. 67, 961 (1991). 

E.G. Drukarev, M.G. Ryskin, and V.A. Sadovnikova, 
larXiv: 1012.03941 



